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We present an exhaustive theoretical analysis of charge and thermoelectric transport in a normal metal- fer¬ 
romagnetic insulator-superconductor (NFIS) junction, and explore the possibility of its use as a sensitive ther¬ 
mometer. We investigated the transfer functions and the intrinsic noise performance for different measurement 
configurations. A common feature of all configurations is that the best temperature noise performance is ob¬ 
tained in the non-linear temperature regime for a structure based on an europium chalcogenide ferromagnetic 
insulator in contact with a superconducting A1 film structure. For an open-circuit configuration, although the 
maximal intrinsic temperature sensitivity can achieve lOnKHz -1 / 2 , a realistic amplifying chain will reduce the 
sensitivity up to 10/rKHz -1 / 2 . To overcome this limitation we propose a measurement scheme in a closed- 
circuit configuration based on state-of-art SQUID detection technology in an inductive setup. In such a case 
we show that temperature noise can be as low as 35nK Hz -1 / 2 . We also discuss a temperature-to-frequency 
converter where the obtained thermo-voltage developed over a Josephson junction operated in the dissipative 
regime is converted into a high-frequency signal. We predict that the structure can generate frequencies up to 
~ 120GHz, and transfer functions up to 200GHz/K at around ~ IK. If operated as electron thermometer, the 
device may provide temperature noise lower than 35nK Hz -1 / 2 thereby being potentially attractive for radiation 
sensing applications. 

PACS numbers: 74.50.+r,85.25.-j,74.25.F-,72.15.Jf 


I. INTRODUCTION 

Recent theories have shown that the spin-splitting induced 
in a superconductor (S) placed in contact with a ferromag¬ 
netic insulator (FI) can be exploited in different kinds of 
spin caloritronic devices such as heat valves 1,2 or thermo¬ 
electric elements 3-6 . They can be used as building blocks 
in phase-coherent thermoelectric transistors 7 , and for the cre¬ 
ation of magnetic fields induced by a temperature gradient 
in Josephson junctions (JJs) due to the thermophase effect 8 . 
NFIS junctions have been also proposed for efficient elec¬ 
tron cooling 9 of the normal metal N 10 . The possible applica¬ 
tions of superconductor-ferromagnetic structures for thermo¬ 
electrics has been also highlighted in a recent review article. 11 

In the present work we theoretically analyze charge and 
thermoelectric transport in a prototype structure based on 
the FIS building block, and explore its application as an 
ultra-sensitive electron thermometer 12-19 and, eventually as a 
temperature-to-frequency converter. Our system consists of a 
normal metal-ferromagnetic insulator-superconductor (NFIS) 
junction, denoted here as the thermoelectric element (TE), 
which is connected, via the superconducting wires Si, to a 
generic load element, as shown in Fig. 1. A temperature 
difference localized between the N and S side of the TE in¬ 
duces a thermoelectric signal 4 . We consider three different 
configurations of the load resistance Rl = °° (open circuit), 
Rl = 0 (closed circuit) and finite load Rl = Rjj where we 
close the system over a generic Josephson element, in the 
dissipative regime, with shunting resistance Rjj. Depend¬ 
ing on the configuration the device will operate in different 
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FIG. 1. (Color online) General scheme of the device based on a nor¬ 
mal metal-ferromagnetic insulator-superconductor (NFIS) junction. 
The latter is shown in the blow-up as stacked layers of different ma¬ 
terials. Ts and 7/y denote the temperature in S and N, respectively, 
Ije is the thermocurrent circulating in the circuit, and V is the ther¬ 
movoltage developed across the thermoelectric element (TE). Si is 
a superconductor contacted to both ends of the TE which contains a 
load with resistance Rl . The latter is intended to be an open-circuit 
Rl —> °°, a dissipationless closed-circuit Rl = 0, or in the form of a 
generic Josephson element operated in the dissipative regime in or¬ 
der to convert the thermo voltage V into radiation at the Josephson 
frequency. 


regimes: i) Seebeck regime , where a Seebeck thermovoltage 
(V) is generated across the TE element at open-circuit; ii) 
Peltier regime , where the gradient of temperature generates 
a circulating thermocurrent that can be probed by an induc¬ 
tive measurement for closed-circuit. Here we explore both 
regimes that includes an estimate of the intrinsic noise and the 
best expected temperature sensitivity with state-of-art technol¬ 
ogy for signal detection. We discuss the advantages and the 
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drawbacks of the different configurations and show that oper¬ 
ated within the non-linear regime, the intrinsic noise of the de¬ 
vice is reduced. In particular, our numerical results show that 
the noise performance is mainly determined by the junction 
differential resistance R 4 which can be drastically reduced be¬ 
yond the linear-response regime with respect to the tempera¬ 
ture. We finally discuss how the generated thermovoltage can 
induce an ac-Josephson effect with a supercurrent oscillating 
at a frequency v = |V|/<f>o 21 , where <f>o ~ 2.067 x 10 -15 Wb 
is the flux quantum. The frequency v can be measured with 
great accuracy providing accurate and fast information about 
temperature difference across the TE. 

The paper is organized as follows: In Sec. II we briefly 
present the general formalism and the expressions for the elec¬ 
tric current flowing through the NFIS junction and the noise as 
a function of all the parameters involved in the system. With 
the help of this expression we analyze in Sec. IIA the elec¬ 
tric and thermoelectric response of the TE in the non-linear 
response regime. In particular, we show the impact of the ex¬ 
change field as well as the role of the barrier polarization on 
the charge current. In Sec. IIB we discuss the different mea¬ 
surement configurations of the device analizing the effect of 
the load resistance Re over the thermo-electrical properties of 
TE recalling the results for the linear regime in Sec. IIC. The 
evaluation of the intrinsic noise properties of the NFIS junc¬ 
tion is done both for the linear an non-linear regime. Assum¬ 
ing a realistic device based on europium sulfide (EuS) as FI 
and superconducting aluminum (Al), operating at low temper¬ 
atures, we discuss the open-circuit and closed-circuit config¬ 
urations respectively in Sec. Ill and Sec. IV. In those sections 
we also discuss the temperature noise performance taking into 
account the most simple measurement scheme with actual 
state-of-art technologies. Finally, in Sec. V, we discuss the 
temperature-to-frequency conversion scheme where the ther¬ 
movoltage developed across the NFIS junction is converted 
into a high-frequency signal by a Josephson element driven 
into the dissipative regime. The full temperature-to-frequency 
conversion capability of the NFIS junction is analyzed, inves¬ 
tigating as well the temperature noise performances. We sum¬ 
marize our results in Sec. VI. 


II. MODEL 


It is instructive to start with the description of the NFIS 
building block. The interaction between the spin of the con¬ 
ducting electrons in S and the localized magnetic moments in 
FI lead to an effective exchange interaction in S that decays 
over the superconducting coherence length £o 26 . We assume 
that the S layer is thinner than £ 0 , so that the exchange field 
(h exc ) induced in S by FI is spatially homogenous. In such a 
case the superconductor density of the states (DoSs) is given 
by the sum of the densities for spin-up (t) and spin-down (|) 
quasiparticles, 
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FIG. 2. (Color online) Characterization of the TE for T§ = T^ (a) 
Current vs voltage {Ije — V) characteristics of the TE element cal¬ 
culated at Ts = 7/v = 0.1 T c , P = 0.5 and for a few values of h exc . 
(b) Differential conductance vs voltage (G — V) characteristics of 
TE calculated for the same parameters as in panel (a), (c) Ije — V 
and (d) G — V characteristics of TE calculated at Ts = 7/v = 0.1 T c , 
h ex c = 0.4Aq and for a few values of P. Aq = 1.764&#T C is the zero- 
temperature, zero-exchange field superconducting gap, T c denotes 
the critical temperature, and Rj is the normal-state resistance of TE. 


Here A(Ts,h exc ) is the pairing potential that depends both 
on temperature Ts in S and h exc , and it is computed self- 
consistently in a standard way 22 from the gap equation 
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where f±(E) = (1 + exp[^+ A 2 =f h exc )\ j , ( 0 D is 
the Debye frequency of the superconductor, Ao is the zero- 
temperature, zero-exchange field superconducting pairing po¬ 
tential, and Jcb is the Boltzmann constant. Furthermore, T 
accounts for broadening, and for an ideal superconductor 

r^o+ 23 . 

We are interested in the current through the NFIS junction 
which in the tunneling limit considered here is given by 4 


Ite = -E jT dE [N+ + PN-] [f N (V,T„) - fs(T s )] . (3) 


Here Rj is the normal-state resistance of the tunneling junc¬ 
tion and N± = (A^ zb N±). Notice that in the tunneling limit 
the Andreev reflection is negligible small and hence no su¬ 
perconducting proximity effect in N takes place. We as¬ 
sume thermalization on both, the S and the N layer neglecting 
any deviation of the distribution functions from their equilib- 
rium form 20 : f s (T s ) = [1 + exp(E/k B T s )]~ 1 and f N (V,T N ) = 
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[1 + exp[(£ + eV)/kBT N ]]~ 1 . Here T N is the temperature in 
the N layer, and — e is the electron charge. The role of the 
FI layer is twofold: it acts as a spin filter with polarization 24 
P = (Gf - G;)/(G^ + G;) and causes the spin-splitting of the 
DoS in the S layer due to the exchange coupling between 
the localized magnetic moments of the FI and the conduct¬ 
ing electrons of S 25-27 . These two features have been demon¬ 
strated in several experiments 28,30-34 . Notice that, according 
to Eq. (3), even in the absence of a voltage bias across the 
junction a finite current Ite can flow provided Tn ^ Ts, as 
demonstrated in Ref. 4. 


A. Electric and thermoelectric response of the TE 

Before analyzing the role of a temperature bias across the 
TE, we determine the current-voltage characteristics (IVCs) 
and differential conductance G = cIIje / dV of the NFIS junc¬ 
tion. We set a low temperature, 7# = Ts = 0.01 T c , where 
T c = Ao/(1.764fcfl) is the critical temperature of the supercon¬ 
ductor. 

The results obtained from Eq. (3) are summarized in Fig. 2. 
Panels (a) and (b) show the I VC and G, respectively, for a 
polarization of the barrier P = 50% and different values of 
the spin-splitting exchange field h exc . In panel (a) one clearly 
sees the deviation of the IVCs from those of a metal-insulator- 
superconductor (NIS) junction. For finite values of h exc there 
is a sizeable subgap current [see Fig. 2(a)] as a consequence 
of the spin-splitting of the DoSs in the S electrode. This 
splitting manifests itself also in the differential conductance 
G [see Fig. 2(b)], where the coherent peaks, usually appear¬ 
ing at V = d=A/ e, are now split in four peaks appearing at 
V = (± A db h exc ) / e. The asymmetry in the height of the coher¬ 
ent peaks stems from the spin polarization P of the FI barrier 
[see Figs. 2(c,d)] where we set h exc = 0.4Ao and the curves 
are calculated for different values of P. Therefore from IVCs 
one can estimate both the polarization of the barrier and the 
spin-splitting induced in S 28 . 

We now assume a finite temperature difference between the 
electrodes 29 , ST = Ts — Tv, and re-calculate the IVCs from 
Eq. (3) for h exc = 0.4Ao and P = 0.9. The results are shown 
in Figs. 3(a) and 3(b), where we keep one of the electrodes 
at temperature 0.0 IT C and change the other electrode temper¬ 
ature. The curves in Fig. 3 reveal two main properties of the 
I VC. First, the I VC strongly depends on the amplitude of the 
temperature difference ST: the larger the temperature differ¬ 
ence, the larger is the current flow at low voltages. In the 
case that the S electrode is heated [see Fig. 3(a)], this trend 
is limited by the reduced critical temperature t; < T c of the 
superconductor originating from the presence of a finite h exc 
which suppresses the A (Ts,h exc ) calculated self-consistently. 
When Ts —» T c *, the TE is driven into the normal state with an 
ohmic characteristic [the red curve in Fig. 3(a)]. 

Second, there is another interesting feature of the IVCs: 
they strongly depend on the sign of ST. For the same value 
of |5T|, the current at V = 0 is larger when the N electrode 
is colder than the S one, i.e., when ST > 0. In other words, 
the thermoelectric effect in the TE strongly depends on the 
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FIG. 3. (Color online) Characterization of the TE for T§ ^ T^ 

(a) Ite — V characteristic of TE calculated for several values of Ts 
(see top legend) at T^ = 0.01 T c , P = 0.9, and h exc = 0.4Aq. (b) The 
same as in panel (a) calculated for several 7/v values at Ts = 0.01 T c . 
Dashed lines in panels (a) and (b) represent the current (///) flowing 
through the Josephson element when it is operated in the resistive 
regime, Ijj = -V/Rjj. 


temperature difference. This feature was not investigated in 
previous works 3,4 in which only the linear response regime 
was discussed. 


B. Measurement configurations 

To use the TE element as a thermometer we need now to ex¬ 
tract from the thermoelectrical signal the temperature gradient 
present across the TE junction. In order to do so we need to 
close the TE circuit over a generic load element modeled with 
a load resistance Re (see Fig. 1). In such a case the voltage V 
developed across the TE for a given ST is the solution of the 
following non-linear integral equation 


Ite (V, T s , T N ,h exc ,P) + 2- = 0, (4) 

Rl 

where Ite is defined in Eq. (3). The solution to the above 
equation is given by the point in which the dashed line (with 
slope proportional to 1 /Rl) in Figs. 3(a) and 3(b) intersects 
the IVCs. 

For a Seebeck-like measurement one needs to maximize the 
thermo voltage opening the circuit, i.e., Rl -a °° and Ite = 0. 
For a Peltier-like measurement, one needs to maximize the 
current closing the circuit with a superconducting loop, i.e., 
Rl = 0 and consequently V = 0. In the case of temperature- 
to-frequency conversion that we will discuss later, one needs 
to include a Josephson element that operates in the dissipa¬ 
tive regime with a load resistance Rl = Rjj , which is the total 
shunting resistance of the Josephson element. 

Independently of the chosen configuration, we assume to 
connect the TE to the detector with two superconducting arms 
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Si. In particular, we assume to place a tunnel barrier between 
S and Si to isolate the S element thereby ensuring its descrip¬ 
tion as a thermally homogeneous superconductor with a spin- 
split DoSs. We neglect here any influence of Si arms such that 
the current through the TE is described by Eq. (3). Supercon¬ 
ductors S and Si are Josephson coupled through the barrier so 
that no additional voltage drop will occur. Furthermore, we 
also assume the NSi junction to be a clean metallic contact, 
thereby contributing negligibly to the total resistance of the 
system and, for simplicity, we disregard the proximity effect 
induced into the N layer by the nearby contacted supercon¬ 
ductor Si 22 . 


C. Linear response regime 


In the linear response regime, i.e., when the voltage V 
and temperature difference ST <C T=(T S + 7]v) /2 across the 
NFIS junction are small, Eq. (4) reads 3,4 . 

8 T 

lTE&lTE = OV+Pa— , (5) 

where 


<7 
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is the electric conductance, and a is the thermoelectric See- 
beck coefficient 4 defined as 


1 

eRj 



EN- 
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(7) 


which, in the linear regime, is connected to the Peltier coef¬ 
ficient n = a/T by Onsager symmetry. Substituting Eq. (5) 
in Eq. (4) and solving respect to the thermovoltage across TE 
element one finds 


yim ^ 


-Pa- 


Rl ST 
R l ° + 1 T 


( 8 ) 


that is valid in the linear response regime assuming a generic 
load resistance Rl- We see immediately that the thermo volt¬ 
age directly measures the temperature gradient in TE. Further¬ 
more, for fixed load resistance, the achievable thermo-voltage 
V increases with the polarization P. In an open-circuit con¬ 
figuration (Rl 00 and Ite = 0) the TE thermo-voltage is 
maximal being 


y7m 


Pa ST 

77 ' 


(9) 


For the closed-circuit (Rl = 0 and V = 0) instead the thermo¬ 
current is maximal being 


I l j n E~Pa—. (10) 

Obviously we see that in the linear regime the open-circuit 
thermovoltage V lm is directly related to the closed-circuit ther¬ 
mocurrent, I lln = oV lin . In particular, the dependence of the 


conversion efficiency on the polarization P and the temper¬ 
ature gradient are the same. This simple picture drastically 
changes if one goes beyond the linear response regime, i.e., 
8 T ~ T . We will see below that the non-linear regime is es¬ 
sential in order to optimize the sensitivity for thermometry 
applications 9 . 


D. Intrinsic noise of TE element 

We now address the zero frequency noise performance of 
the NFIS junction. In this case, the main source of noise is 
the current noise (Si) generated in the TE that is described by 
generalizing the expression derived in Ref. 37 in the presence 
of a ferromagnetic tunneling barrier: 

s, = E J~JE[N + +PN-]M(E,V,T n ,Ts), (11) 


where 

M = f N (VJ N )[l-fs(T s )] +fs(T s )[l-MV,T N )i (12) 

and the bias V is given by the solution of Eq. (4). We note that 
the previous formula describes both thermal, i.e., eV <C Tv, Ts, 
and shot noise, i.e., T/v,7^ <C eV, and holds in the tunneling 
regime. 

Previous expression simplifies in the linear response regime 
discussed before where we can neglect any term 0(ST) in 
Eq. (11) finding the thermal noise 

S l / n = E \N+ +PN-]f N (T)[l -f s (T )}, (13) 

which may be expressed as 

S l / n =4k B Ta, (14) 

where o is the TE electric conductance of Eq. (6). In the open- 
circuit configuration, it is more convenient to write the voltage 
noise spectral density 

S l j n = 4 Jc b T/o . (15) 

Below we show that in the non-linear regime one can approx¬ 
imate Eqs. (14)-(15) by substituting o by l/Rd, where Rd is 
the TE differential resistance. 


III. TEMPERATURE-TO-VOLTAGE CONVERSION 

When TE is in an open-circuit configuration (Rl °°), one 
can realize a temperature-to-voltage conversion scheme. In 
such case no charge current flows through the TE, 

I T E(Vo,Ts,TN,h exc ,P)=0 , (16) 

then a voltage Vo develops across the TE for ST / 0. The 
value of Vo can be obtained from the solution of Eq. (16). The 
results are shown in the two upper panels of Fig. 4. Specif¬ 
ically, panel 4(a) shows the dependence of Vo on Ts for dif¬ 
ferent values of h exc , P = 0.9 and 7# = 0.01T C . The increase 
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FIG. 4. (Color online) Open-circuit configuration, (a) Thermo- 

voltage Vo vs Ts calculated for a few values of h exc (see legend) at 
7/v = 0.01 and P = 0.9. (b) Vo vs Tn calculated for a few values of 
h exc at 7^ = 0.01T C and P = 0.9 (full lines), or P = 1 (dashed lines), 
(c) and (d) show the absolute value of the corresponding temperature- 
to-voltage transfer function \xy\ calculated for the same values of 
hexc • 


FIG. 5. (Color online) Open-circuit configuration, (a) Thermo¬ 
voltage Vo vs T s calculated for a few values of T (see legend) at 
T n = 0.01 T c , h exc = 0.3Ao, and P = 0.9. (b) Vo vs Tn calculated 
for a few values of T at Ts = 0.01 T c , h exc = 0.3Ao, and P = 0.9 
(full lines), (c) and (d) show the absolute value of the corresponding 
temperature-to-voltage transfer function | t y | calculated for the same 
values of T. 


of Ts, from the minimal temperature T N , leads first to an en¬ 
hancement of | Vo | - A further increase of Ts leads to the sup¬ 
pression of the superconducting energy gap and a correspond¬ 
ing suppression of Vo- The voltage Vo vanishes when when su¬ 
perconductivity is fully suppressed for Ts T* . We note that 
Vo reaches zero continuously owing to the fact that we have 
chosen values of h exc for which the superconducting-normal 
state transition is of the second order 38 

A different temperature behavior of Vo is obtained when S 
is kept at Ts = 0.01 T c and Tn is varied, as shown in Fig. 4(b). 
In particular, besides the obvious change of sign, Vo grows 
monotonically by increasing Tn until it reaches an asymptotic 
value. It is important to stress that the curves Vo(7/v) depends 
strongly on the polarization P of the barrier [see Fig. 4(b)]. In 
particular, the larger P the larger is the thermovoltage Vo(Tn) 
developed across the TE. By contrast, the Vo (7^) amplitude 
turns out to be almost unaffected by the value of P. 

The different behaviors as a function of 8 T allow one to 
reconstruct both the amplitude and direction of the thermal 
gradient in the TE element. This further information could 
be eventually exploited to reconstruct the spatial position of 
a heating event, thereby opening interesting possibilities to 
build detector-like devices. 

An useful figure of merit to estimate the performance of 
the TE is the temperature-to-voltage transfer function, Ty = 


dV/dT. The absolute value of this quantity is shown in 
Fig. 4(c) for ST > 0 and in Fig. 4(d) for ST < 0. We have 
normalized it to the natural unit Ao/ eT c . In Fig. 4(d) we show 
the case of barrier polarization selectivity P = 1 which corre¬ 
sponds to the case with maximal possible transfer function at 
given h exc . 

In order to show the impact of the broadening parameter we 
display in Fig. 5 the same quantities as in Fig. 4 but calculated 
for fixed h exc = 0.3Ao, P = 0.9 and for different values of T 
ranging from 10 -6 Ao up to 10 -2 Ao 39-41 . The overall qualita¬ 
tive behavior and the order of magnitude of the effect is the 
same for all these values. From a quantitative point of view, 
the temperature-to-voltage conversion turns out to be less ef¬ 
fective the larger the value of T. Throughout the paper we as¬ 
sume r = 10 _4 Ao which is the typical value for conventional 
Al-based superconducting junctions 9,39 . 


A. Noise performance analysis for the open-circuit 
configuration 

We now focus our analysis on the noise performance of 
the temperature-to-voltage conversion with the NFIS junction. 
We need to convert it in a voltage noise assuming that load re¬ 
sistance Rl 00 . This means that the voltage noise spectral 
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FIG. 6. (Color online) Open-circuit configuration, (a) Square 

1/2 

root of the voltage noise spectral density Sy vs T$ calculated at 

Tn = 0.0 ir c for a few values of h exc (see legend), (b) Sy 2 vs 7# 
calculated at T$ = 0.01 T c for the same h exc values, (c) Temperature 
sensitivity st vs T$ calculated for the case of panel (a), (d) sj vs 7# 
calculated for the case of panel (b). Dashed lines show st as calcu¬ 
lated using the approximate formula Eq.19 for h exc / Aq = 0.5. In all 
these calculations we set P = 0.98, and T c = 3K. 


density (Sy) generated across the TE, is 36 

S v (V 0 .Ts.T N ..h eX c,P) = S,R 2 d .. (17) 

where the R d = dVo/dljE is the differential resistance of the 
TE, and the bias Vo is given by the solution of Eq. (16). 

In Fig. 6(a)-(b) the square root of noise spectral density 
i^/Sy) is displayed for a TE element with a barrier character¬ 
ized by a realistic value of polarization P = 0.98. This spin- 
filter efficiency is representative for EuO or EuS FI barriers, 35 
and we assume a superconductor with T c = 3K which would 
be implementable with ultra-thin A1 films 30-33 . 

We find that for 8 T > 0 the minimal noise value is obtained 
in the non-linear regime 7/y <C Ts < T c where the voltage noise 
can be as low as ~ 600fV Hz -1 / 2 , and is two orders of mag¬ 
nitude lower than the equilibrium case 8 T/T <C 1, where 
T = (Ts + Tn )/2 is the average temperature. For 8 T < 0, the 
noise performance is worse being at best a few tens of pV 
Hz -1 / 2 for the non-linear regime Ts ^ T c . 

The intrinsic temperature noise (temperature sensitivity) 
per unit bandwidth of the thermometer (sj) is related to the 
voltage noise spectral density as 


st = 




M' 


(18) 


In Fig. 6(c)-(d) we show the temperature noise sj for the 
open-circuit configuration for the two cases of Fig. 6(a)- (b). 
The differences between the voltage spectral density is en¬ 
tirely given by the transfer function which is highly non-linear 
as a function of 8 T. We notice that in the linear regime 
\8T\/T < 1 the temperature noise, given by Eq. (15), is a 
few tens of /iK Hz -1 / 2 . The maximum temperature sensibil¬ 
ity is obtained in the nonlinear regime, i.e., 18 T\/T 1 where 

the temperature noise can be as low as 8nK Hz -1 / 2 coincid¬ 
ing with the minimal voltage noise. By contrast, for the case 
8 T < 0, the best noise performance is around 180nK Hz x / 2 . 

It is interesting to observe the scaling behaviour of noise 
power as function of the junction normal-state resistance Rj. 
Indeed as Si oc \/R T [see Eq. (11)], from Eq. (17) one can 
conclude that Sy oc R T since R d °z Rt- At the same time there 
is no scaling behaviour of the transfer function, t y = dVo/dT. 
This may be easily inferred, for instance, from the expression 
of the thermovoltage Vo in the linear regime, Eq. (9), since 
V lin oc a/a which is the ratio of two quantities with the same 
scaling \/Rj. Alternatively one can deduce it from the rela¬ 
tion between open-circuit voltage an the temperature differ¬ 
ence, 8 T Eq. (16), where Rj enters only as an overall pref¬ 
actor. We conclude that sj y/Rr which shows immediately 
that the reduction of TE resistance would be beneficial for in¬ 
creasing the sensitivity in temperature measurement. 

These considerations suggest that in the non-linear regime 
the differential resistance R d takes the role of Rj. Indeed one 
can guess a way to generalize Eq. (15) to the non-linear regime 
by replacing the linear conductance a by the differential con¬ 
ductance l/Rd such that 

Sy ~ 4k B TR d (19) 

where the temperature is taken as the average T = (Ts + 
7]v)/2. The previous expression should converge to the lin¬ 
ear result when \8 T\/T <C 1. The full numerical results of 
Fig. 6(c)-(d) demonstrate the accuracy of Eq. (19) shown as 
dashed lines. This shows that the noise performance is es¬ 
sentially characterized by the dependence of differential re¬ 
sistance R d . Therefore this approximation is extremely useful 
to estimate the noise performance with the knowledge of the 
differential resistance R d only. 

It is important to emphasize that in a realistic measurement 
scheme the temperature sensitivity of the device is limited 
by the amplifying chain. Indeed, in general, the voltage sig¬ 
nal must be amplified with a low-noise preamplifier which is 
characterized by its intrinsic voltage noise. The preamplifier 
noise may degrade the total noise performance. In particu¬ 
lar, assuming for the pre-amp a square root spectral density 
of ^ InV Hz -1 / 2 , it is clear that it will dominate over the in¬ 
trinsic voltage noise of the signal which can be smaller by a 
few orders of magnitude [see Fig. 6(a)-(b)]. The preamplifier 
is therefore the main bottleneck to the temperature detection 
in this configuration scheme although it has the advantage of 
suppressing the noise non-linearities over the considered tem¬ 
perature window. The realistic performance of this measure¬ 
ment scheme will be roughly ~ 10/lK Hz -1 / 2 . This limitation 
can be overcome by exploiting a close-circuit configuration, 
as discussed in the next section. 
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FIG. 7. (Color online) Closed-circuit configuration (a) Thermocur¬ 
rent Ij E vs Ts calculated for a few values of h exc at 7/v = 0.01 T c . (b) 
I EE vs 7/v calculated for the same values of h exc . (c) and (d) show the 
absolute value of the temperature-to-current transfer function |t/| vs 
temperature for the panel (a) and (b), respectively, calculated for the 
same values of h exc . In all these calculations we set P = 0.98. 
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FIG. 8. (Color online) Closed-circuit configuration (a) Square root 

1 /2 

of the current noise spectral density S/ vs Ts for 7/v = 0.01 T c and 

1/2 

for different values of h exc (see legend), (b) S/ vs 7/v for Ts = 
0.0 IT C for the same values of h exc . (c) st vs Ts for the case of panel 
(a), (d) st vs 7/v for the case of panel (b). Dashed lines shows sj as 
calculated using the approximate formula Eq. (22) for h exc / Aq = 0.5. 
In all these calculations we set P = 0.98 and T c = 3K. 


IV. TEMPERATURE-TO-CURRENT CONVERSION 

Hereafter we analyse the performance of a closed-circuit 
configuration which correspond to temperature-to-current 
conversion. In this setup the TE current Ije (V = 0) = I EE 
is given by Eq. (3) which depends only on Ts and T N . 

In Fig. 7(a) we show how the current depends on Ts for 
different values of h exc keeping fixed the barrier polarization 
P = 0.98 and 7# = 0.01T C . The general behaviour has ’’shark- 
fin” shape which increases in amplitude with the h exc . After 
reaching a maximum at T£ the I EE decreases with Ts until the 
critical temperature T* is reached and the superconductivity 
is completely suppressed. 

If we fix Ts = 0.017A by changing 7# we get for 8 T < 0 an 
obvious opposite sign for the thermocurrent I EE , and an ab¬ 
solute value of the thermocurrent \I EE \ which monotonously 
increases by enhancing |5T|. It finally saturates to the maxi¬ 
mal value 

ITE^ax = 2^ f_J E \- N + +PN ~} S S n ( £ )’ ^ 20) 

which is easily obtained from the general expression of the 
TE current, Eq. (3), by taking the limit 7^ -> 0 and T N 
oo with V = 0. From this, we can conclude that an arbitrary 


enhancement of |5T| is not of particular benefit to increase 
the current signal. 

In Fig. 7(c)-(d) we show the absolute value of the 
temperature-to-current transfer function T/ = dI EE /dT re¬ 
spectively for the case (a) and (b) of the same figure. For 
8 T > 0, the transfer function has two different behaviors de¬ 
pending if Ts is smaller or larger than T £. On the other hand 
one sees that, independently of the sign of ST, the transfer 
function is maximized in the non-linear regime |5r| ~ T. 
This is an important issue in order to increase the tempera¬ 
ture sensitivity. 


A. Noise performance analysis for the closed-circuit 
configuration 

In Fig. 8(a)-(b) the current noise 5/ of the closed-circuit 
configuration is shown as obtained from Eq. (11) with V = 0. 
The current noise as a function of ST is minimized in the lin¬ 
ear regime obtaining ~ 15fAHz -1 / 2 , and grows by increasing 
|5r|. The noise behaviour of the closed-circuit configuration 
it less affected by the sign of 8 T in comparison to the open- 
circuit one (see Sec. Ill A). The current noise increases with 
h exc since also the average current I EE increases in such a case 
[see Fig. 7(a) and (b)]. 
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The intrinsic temperature noise per unit bandwidth of the 
thermometer st in this configuration is given by 



( 21 ) 


where |tj| is the temperature-to-current transfer function dis¬ 
cussed before. The same scaling behaviour shown before for 
st y/Rr still holds in this configuration since now 5/ oc 1 /R T 
but T/ oc 1 /R t . Consequently also for this case the minimiza¬ 
tion of Rt would be, in general, beneficial for improving noise 
performance. As in the previous section one can try to gen¬ 
eralise this argument for the non-linear regime by replacing 
Rt with Rj. Since Rj is largely reduced in comparison to the 
linear regime value l/a, one can expect an increase of the 
noise in the non-linear regime. At a first glance, this does 
not look plausible since the current noise is in general higher 
[see Fig. 8(a)-(b)]. However, our guess seems to be correct, 
as shown in Fig. 8(c) and (d), where the temperature noise is 
minimized for large values of 8 T. 

The lowest intrinsic noise ^ 35nK Hz -1 / 2 is obtained in the 
non-linear regime for 8 T « IK, when P = 0.98 and T c = 3K 
are chosen. The main difference with the open-circuit config¬ 
uration (see Fig. 6) is that in the present situation the noise 
depends weakly on the sign of 8 T in a wide temperature re¬ 
gion. Moreover, in contrast to the open-circuit configuration, 
the noise shows a rather smooth behavior. 

As pointed out above, the behaviour of the current noise in 
the non-linear regime can be approximated by the expression 


Si 


W 

Rrl ’ 


( 22 ) 


where the linear conductance a of Eq. (14) is replaced by 
the inverse differential resistance, l/Rd- In Fig. 8(c) and (d) 
we show (dashed lines) this approximation for the case where 
we expect the largest nonlinearities, i.e., for h exc /Ao = 0.5. 
We can thus conclude that this simple formula gives a fairly 
accurate description of Si in the non-linear regime. 

In terms of overall temperature noise the closed-circuit con¬ 
figuration has two advantages: First, the smooth behaviour 
of temperature noise makes it more attractive than the open- 
circuit configuration. Second, while the ideal noise is better 
for the open-circuit configuration (see Figs. 6 and 8), one 
needs to evaluate the total noise of the measurement which 
includes the addition of the preamplifier noise. The latter, as 
discussed in the previous section, strongly degrades the result¬ 
ing noise figure. The close-circuit configuration offers a way 
to overcome this limitation, as discussed in the following. 

Specifically, we propose to measure the current signal by 
coupling the closed circuit via a mutual inductance M to a 
SQUID. The latter measures the flux generated by the current 
circulating in the thermoelectric circuit. The total tempera¬ 
ture sensitivity, which includes now the SQUID noise, can be 
written as 


_ v ,s '+( x T"> 2 > 

% W It,I 


( 23 ) 



FIG. 9. (Color online) Characterization of the temperature-to- 
frequency converter (a) Frequency v vs Ts calculated for a few val¬ 
ues of h exc at 7/v = 0.017^. (b) v vs at Ts = 0.01 T c calculated 
for the same h exc values as in panel (a), (c) and (d) show the ab¬ 
solute value transfer function |t v | of panel (a) and (b), respectively, 
calculated for the same values of h exc . In all these calculations we 
setP = 0.98, Rt/Rjj = 0.2, and T c = 3K. 


where the temperature-to-flux transfer function is = Mi/, 
and the TE flux spectral density S^ E = M 2 Si is added to the 
SQUID noise (S S ® UID ) to give the total flux noise, S^ OT = 
Sq E +S S ® UID . The square root of the flux noise for high- 
quality commercial SQUID can be as low as yjS S ® UID ~ 

10 -7 <f>o Hz -1 / 2 , which is then converted into an effective cir¬ 
culating current noise in the thermoelectric circuit of ~ 20 
fA Hz -1 / 2 by dividing it with a typical value for the mutual 
inductance M = 10 _8 H. By looking at Fig. 8(a) and (b) we 
immediately see that the intrinsic TE current noise will, in 
general, dominate over the SQUID noise almost everywhere 
in the non-linear regime where we can achieve the best sensi¬ 
tivity. Therefore, the temperature noise of this measurement 
scheme is only limited by the intrinsic TE noise mechanisms, 
and can be as low as ^ 35nK Hz -1 / 2 for a moderate tempera¬ 
ture non-linearity (see Fig. 8). 


V. TEMPERATURE-TO-FREQUENCY CONVERSION 

We now focus on the temperature-to-frequency conversion 
process. This conversion is achieved with the device sketched 
in Fig. 1 where the thermovoltage generated across the TE is 
applied to a generic Josephson element. The latter is set to op- 
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erate in the dissipative regime when Ijj = V/Rjj, where Rjj is 
the total shunting resistance of the Josephson element. In this 
case, there is a time-oscillating current through the Josephson 
element with a frequency equal to the Josephson frequency, 
v = |V|/O> 0 - As discussed above, the value of V depends on 
the temperature difference 57 across the TE, and therefore the 
frequency emitted by the Josephson junction is a measure of 
ST. 

In order to quantify the temperature-to-frequency conver¬ 
sion effect, one has to determine the voltage V developed for 
any given ST imposed across the TE which satisfies Eq. (4) 
with a finite load resistance Rl = Rjj- This configuration is 
intermediate between the open-circuit and the closed-circuit 
setup discussed above. In the following, we consider the 
case where Rj /Rjj = 0.2 in order to produce a detectable fre¬ 
quency signal between 10 GHz and fractions of THz. 

The frequency-to-temperature performance of this configu¬ 
ration is shown in Fig. 9. We again used the spin-filter effi¬ 
ciency P = 0.98 and the critical temperature T c = 3K adopted 
in the previous sections. Panels 9(a) and 9(b) show the fre¬ 
quency generated by the Josephson element, for positive and 
negative ST, respectively. In the linear response regime, the 
TE thermovoltage depends only on |57| as can be seen from 
Fig. 9. The information about the sign of ST is eventually 
recovered only for the non-linear regime. 

If 7yv is kept at 0.01 T c the maximum frequency is achieved 
around h exc « 0.2Ao for Ts « 0.757 c , and obtains values as 
large as ~ 120GHz. If Ts is kept at low temperature, v in¬ 
creases monotonically by increasing both 7# and/or h exc and 
obtains a maximum of ~ 80GHz 

In the present setup an important figure of merit of the struc¬ 
ture is represented by the temperature-to-frequency transfer 
function, T v = dv/dT , plotted in absolute value in Figs. 9(c) 
and 9(d). In particular, |t v | exceeding 200GHz/K around 
Ts ~ IK can be achieved for h exc = 0.5Aq by heating S, while 
|t v | up to ~ 55 GHz/K can be achieved with the same values 
by heating N. 


A. Noise performance 


In the temperature-to-frequency conversion process, the 
noise is determined by the bias fluctuations generated from the 
current noise via the load resistance seen by the TE, i.e., the 
parallel between the Josephson element total resistance Rjj 
and the TE resistance Rd’. Tl = RdRjj/(Rd + Rjj)- Note that 
the differential resistance Rd = dVo/dljE is calculated from 
the solutions of Eq. (4) where Rl = Rjj- The important quan¬ 
tity is represented by the frequency noise spectral density (5 V ) 
which can be expressed as 


Sy — 


Sj'R 2 


(24) 


Finally, the intrinsic temperature noise per unit bandwidth of 
the thermometer (st) is related to the frequency noise spectral 
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FIG. 10. (Color online) Noise performance of the temperature-to- 
frequency converter (a) Square root of the frequency noise spectral 

1 /2 

density Sy vs Ts calculated at 7/v = 0.017 C for a few values of h exc . 

1/2 

(b) Sy vs 7/v calculated at Ts = 0.01 T c for the same h exc values. 

(c) Temperature sensitivity sj vs Ts calculated at 7/v = 0.01 T c for a 
few values of h exc . (d) sj vs 7Jv calculated at 7^ = 0.017 c for the 
same h exc values. In all these calculations we set P = 0.98, T c = 3K, 
Rt/Rjj = 0.2, and Rj = IQ,. 


density as 


s T = 



(25) 


Figure 10(a) and (b) show the calculated square root of 
the frequency noise spectral density S v for positive and neg¬ 
ative 57, respectively, calculated for the same parameters as 
in Fig. 9, and for Rt = 1£2. In particular, for positive 57, 
the noise spectrum S v shows a non-monotonic behavior with 
a maximum at intermediate temperatures, and suppression at 
higher 57. By contrast, for 57 <0 the noise spectrum grows 
monotonically with 1571, and it is less influenced by h exc . 

The behavior of st is displayed in Fig. 10(c) and (d). At 
small 1571, in the linear regime, the noise sensitivity is given 
by several tens of /iK Hz -1 / 2 . By increasing 1571 the growth 

of Sy [Figs. 10(a) and (b)] is advantageously compensated 
by the enhancement of T v [see Fig. 9(c) and (d)]. The best 
noise performance is obtained when |t v | it is quite near its 
maximum. The values of st ~ 35 nK Hz -1 / 2 is obtained 
around IK for h exc = 0.5Aq. After the minimum of st , for 
57 > 0, we see a peak due to the divergence of |t v | _1 , i.e., 
the vanishing of the transfer function, as shown in Fig. 9(c). 
Differently, for 57 <0 one observes a smooth increase of st 
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determined by the progressive reduction of the transfer func¬ 
tion [see Fig. 9(d)] which is consequence of the saturation of 
the frequency when 7# °°. We conclude by noticing that, 

also for this case, the best noise performance is obtained in 
the non-linear regime for |<57| < IK where St is almost inde¬ 
pendent of the sign of 8 T. This is essentially due again to the 
fact that Rd is strongly reduced in the non-linear regime. 

In this configuration the power of the generated frequency 
signal might be somewhat low. Anyway, one can deploy the 
standard techniques in order to increase the emission power 
by connecting in parallel arrays of JJs. 42,43 . The only limit¬ 
ing factor in that case would be the power that the TE element 
could sustain and transfer to the JJs. A rough estimate shows 
that when Rt/Rjj = 0.2 and Rt = 1 kl the TE could produce 
a power of the order of ~ 100pW...10nW which would be 
high enough to make the 10 — 100GHz signal generated by 
the Josephson junction detectable. 

VI. SUMMARY 

In summary, we have theoretically investigated a ther¬ 
moelectric structure based on a normal metal-ferromagnetic 
insulator-superconductor (NFIS) junction. We fully charac¬ 
terize the thermoelectrical properties of the TE both in the 
linear and non-linear regimes. We assumed different mea¬ 
surement configuration as determined by the load resistance 
value. In particular, we showed that by exploiting realistic 
materials such as EuS or EuO (providing polarization P up to 
~ 98%) in combination with superconducting A1 thin films 
the device is able to provide remarkable temperature noise 
performances. We find that in the open circuit configura¬ 
tion, where the temperature signal is returned via the See- 
beck thermo-voltage, the lowest achievable intrinsic noise of 
~ lOnKHz -1 / 2 is limited by the amplifying chain. On the 
other side, we found that in the closed-circuit configuration, 
where the temperature information is encoded in the Peltier 
thermo-current , one can detect the signal via a low-noise flux 
measurement of an inductively-coupled SQUID. In such case 
the temperature noise performance is mainly determined by 
intrinsic noise mechanisms, with the best value of ~ 35nK 
Hz -1 / 2 achievable with state-of-art SQUID technology. In¬ 
terestingly, we identified in the differential resistance Rd of 
the TE one of the main factor that determines the intrinsic 


noise performance of the system. This explain why the best 
noise performances are obtained in the non-linear temperature 
regime since for that regime Rd is strongly suppressed. This 
is a non-trivial consequence of the strong non-linearities pe¬ 
culiar of the NIFS junction. 

We finally discuss a temperature-to-frequency converter 
where the obtained thermovoltage is converted through a dis¬ 
sipative Josephson junction into a high frequency signal in the 
frequency window spanning from a few GHz up to 10 u Hz. 
In particular, we have shown that the device allows for the 
generation of Josephson radiation at a frequency that depends 
on both the amplitude and sign of the temperature difference 
across the NFIS junction therefore opening the route for high- 
frequency detection associated to high temperature sensitivity. 
Frequencies up to ~ 120GHz and large transfer functions (i.e., 
up to 200GHz/K) around ~ 1 — 2K can be obtained in a struc¬ 
ture implementable with the above mentioned prototype FIs. 
In this configuration the device is capable to provide intrinsic 
temperature noise down to ~ 35nKHz -1 / 2 around IK for a 
sufficiently large h exc . The proposed superconducting hybrid 
structure has the potential for the realization of effective on- 
demand on-chip temperature-to-frequency converters as well 
as ultrasensitive electron thermometers or radiation sensors 
easily integrable with current superconducting electronics. 
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